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We employ Chiral Perturbation Theory (ChPT) to evaluate the complete set of pseudoscalar 
matrix elements, < 0j_P/ |7r° , 77 >, with Pj any of the flavor-diagonal pseudoscalar currents (/ = 
u, d, s), to 0(md — m u ), and to next-to-leading order in the chiral expansion. These matrix elements 
represent the basic input to a QCD sum rule analysis of isospin breaking in the nNN couplings 
using the three-point-function method. We discuss also how one could use the results to construct a 
one parameter family of interpolating fields for the 7r°, all of whose members have zero vacuum-to- n 
matrix element, and explain how this could in principle be used to provide non-trivial tests of the 
reliability of the assumptions underlying the use of the three-point-function method. It is shown that 
the isospin-breaking mixing parameters required for this construction receive significant corrections 
beyond leading order in the chiral expansion. 
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I. INTRODUCTION 



QCD sum rules [1] treatments of the isospin-conserving nNN coupling, both using the two-point-function and three- 
point-function methods [2-8], have proven quite successful. As a consequence, it is natural to consider extending these 
treatments to an investigation of the strong interaction (m^ — m u ^ 0) contributions to the isospin breaking of these 
same couplings. The resulting differences among g pp7T o , g nn7r o and <? np7r + would provide contributions to isospin- 
breaking observables in few-body systems, and an evaluation based on methods closely related to QCD is, therefore, 
highly desirable, particularly in light of the recent revival of interest in such phenomena (especially charge symmetry 
breaking (CSB)) [9-25]. 

A first attempt in this direction was made in Ref. [25] (MH), where a sum rule analysis using the three-point-function 
method was performed. In this approach, one considers the three-point functions 

A N >Nn(pi,P2,q) = J d i xid 4 x 2 exp[i(pi.xi - p 2 .x 2 )} < 0\T (r) N i (xi)P„(0)f} N (x 2 )) |0 > (1) 

where rj^^' are (Ioffe) interpolating fields for the nucleons N', N and P^ a pseudoscalar interpolating field for the pion 
of interest (which carries momentum q — pi — p 2 ). The desired coupling, gN'N-n, then appears as a coefficient in the 
contribution to the phenomenological ("hadronic") representation of A^'Ntt having simultaneous poles at p\ = p\ = 
Mjy and q 2 = m 2 . To evaluate these couplings in terms of vacuum condensates, one then simultaneously evaluates 
An'Nk using the operator product expansion (OPE) (the validity of which requires q 2 < —1 GeV 2 ), and relates the 
hadronic and OPE representations of the correlator via an appropriately Borel transformed dispersion relation [1]. 
The contributions due to the pion coupling are isolated by keeping only those terms in the OPE proportional to ^75 
(as in the phenomenological representation) which have a 1 / q 2 pole in the limit that the quark masses arc neglected 
[2-4,25]. 

It is necessary to exercise some care in the choice of the interpolating field for the neutral pion in Eq. (1). Indeed, 
if one were to make the most obvious choice, i.e., the 1=1 combination 

= Pi=i = P U -Pd (2) 

where Pf = fvysf is the pseudoscalar current for quark flavor, /, one would find that because, to 0{rrid — m u ), P n3 
also couples to the 77, the analysis for the extraction of the isospin-breaking ttNN couplings was contaminated by 77 
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contributions. This contamination is of the same order in isospin breaking as the ttNN isospin-breaking couplings of 
interest, and hence not negligible [25]. This problem was overcome in MH by the alternate choice 

P^ H = (Pu-P d ) + e MH P 8 , (3) 

where 

Ps = -j= (Pu + Pd- 2P S ) (4) 

and e MH is to be chosen such that < 0\P£o H \r) >= 0. To leading order in the chiral expansion, one has, from the 
results of Gasser and Leutwyler [26], 

MH a V3 f m d -m u 

£ = Vci = I 

where 9q is the leading order 7r° — rj mixing angle. Since, however, it is known that next-to-leading order contributions 
to, e.g., f n are large (~ 30%), it is likely that similar next-to-leading order corrections to e MH will also be important. 

In this paper we will compute the pseudoscalar current matrix elements, < 0|P/|7T°,7? >, required to evaluate 
these corrections, to 0(md — m u ) and to ncxt-to-leading (1-loop) order in the chiral expansion. We will, in addition, 
show that it is possible to use these matrix elements to construct, not just the field P^i H above, but a continuous 
family of similar interpolating fields, all with zero vacuum-to-r; matrix element. These fields can then, in principle, 
be used to study the reliability of the assumptions which go into the use of the three-point-function treatment of the 
isospin-breaking itNN couplings, a topic we will return to in a later publication. 

Let us consider, now, the problem of constructing "pure" ir° and rj interpolating fields from the flavor diagonal 
pseudoscalar currents. If P is any 1 = pseudoscalar field combination having non-zero 77 matrix element, then 
forming 



t MH = g = - (5) 

4 V m s — m ' 



and requiring < dlP^o^ >= 0, i.e., 



P„o = (P u - P d ) + e n P (6) 
. < 0\P U - Pd\v > 



<0\P \V> 



(7) 



produces an interpolating field, P^o, which couples to the physical 7T°, but not at all to the r\. A "pure" 77 interpolating 
field is constructed analogously: 

Pr, = Po + er,(Pu ~ Pd) (8) 

with 

<0\P u -P d \*°>- W 

The basic ingredients in the construction of P n o, P v are, therefore, the vacuum-to- 7r° , r\ matrix elements of the 1=1 
(Pr— 1) and 1 = (P u +Pd and P s ) pseudoscalar currents, where 7r°, 77 refer to the physical (isospin-mixcd) pseudoscalar 
mesons. It is these matrix elements that we will evaluate to 0(md — m u ), and to 1-loop order in Chiral Perturbation 
Theory (ChPT) (next-to-leading order in the chiral expansion) below. 



II. EVALUATION OF THE VACUUM-TO- tt°, r/ MATRIX ELEMENTS OF THE PSEUDOSCALAR 

CURRENTS 

Having explained above why an evaluation of the isospin-breaking pseudoscalar current vacuum-to- 7r° , 77 matrix 
elements is of interest, we now turn to the explicit evaluation of these quantities, for which we employ the methods 
of ChPT [26]. In order to obtain both the leading and next-to- leading contributions in the chiral expansion of these 
quantities, it is necessary to include those pieces of the effective Lagrangian up to and including fourth order in the 
chiral counting [27], in the presence of scalar and pseudoscalar external sources, s and p. The relevant terms, as 
worked out by Gasser and Leutwyler [26], are 
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£ eff = £^ + C& 



(10) 



where 



and 



C {A) = L x [Tiid^Ud"^)] 2 + L 2 r Tt{d lx Ud v U^) r It{d^Ud v U^) + L 3 Tr(<9 M C/ 1 d"Ud u U^ d v U) 
+ L 4 Tr(d^Ud^)Tr( X f U + x*7 f ) + L 5 Tr[d^d^U( X f U + U^ X )] 
+ U [Tr&U + xU^)} 2 + L 7 [Tt( X ^U - X U^)] 2 + L 8 Tt(^U X ^U + x^X^) 
+ H 2 Tr(xh) 



(11) 



(12) 



where U = exp(iA.7r//), with {7r a } the octet of pseudoscalar Goldstone boson fields and {A a } the usual Gell-Mann 
matrices, / is a (second order) low-energy constant (LEC) which turns out to be equal to the pion decay constant in 
the chiral limit, 



X = 2B (s + ip) 



(13) 



with B another second order LEC, related to the quark condensate in the chiral limit, and {L k }, H 2 are fourth 
order LEC's. Terms in the general fourth order piece of the effective Lagrangian, £ e ff, which vanish in the absence 
of external vector and axial vector sources have been dropped in writing Eqs. (11), (12). The second and fourth order 
LEC's appearing in these equations are not determined by the symmetry arguments which go into constructing jC e ff 
and must be determined phenomenologically. (See, for example, Ref. [28] for a recent up-to-date tabulation). For our 
purposes, the external scalar source is to be set equal to the quark mass matrix, M, in order to correctly incorporate 
the effects of explicit chiral symmetry breaking, and the external pseudoscalar source is to be retained to first order, in 
order to generate the correct representation of the underlying QCD pseudoscalar currents in the low-energy effective 
theory. 

From the form of £ e // given above, it is straightforward to obtain the following low-energy representations of the 
pseudoscalar currents in terms of the Goldston boson fields {n a } occuring in U : 



Pu = B f 



7T3 + 



7T8 

V3 




3^8 



2 

7T37r 8 



+ 27T+^ 3 

3v3 



-—tt+tt-tts + 4K+K-n 3 + 4K K°TT 3 + 2V2n + K-K° + 2V2tt-K + K° 



+ 



32B% 



I 

+m u I tt 3 



7T8 



V3 



(m s + 2m) I 7T 3 + —= ) L 6 + (m u - m d )-K 3 + —={m- m s )-K S L 7 



V3 V 



TTs_ 

V3 



L, 



+ 



Pd = B [ 



3 7r 8 - 7T 3 7T 8 + 



3V3 



27T+7T" 



7T3 



+ Tr+Tr-Trg _ 4K+K-n 3 - 4K°K°t: 3 + 2V2ir+K-K° + 2V2V K+K° 
V3 



+ 



32.Bg 



^8 



V3 



(m s + 2m) -7T3 + —j= ) L 6 + [ (m u - m d )-n: 3 + —j={m - m s )7r 8 L 7 



V3 V 



+m d -7T 3 + 



P s = B f 




2K+K~ tt 3 - 2K K a n 3 - 2K°K°tt 3 - -^K+K'tts 

v3 



3 



-^=K°K°n 8 + 2V2ir+K-K° + 2y/2n~ K+ K° 



+ 



32B 2 



-—7=(m s + 2m)Tr a L 6 
v3 



+ (m u - m d )ir 3 + — = (m - m s )n s L 7 -=m s -K%L 



73 



(14) 



where ir 3 and irg are the pure 1—1,0 unrenormalized octet fields appearing in the matrix variable U , rh = (m u +md)/2, 
and we have written down explicitly only those terms which will be required for our calculations below. 

Using these representations, it is then straightforward to compute the matrix elements in question. The set of graphs 
which contribute is shown in Fig. 1. For practical purposes it is convenient to recast the low-energy representation 
of the currents in terms of the renormalized, physical ir°, r\ fields, rather than unrenormalized octet fields tt 3 and tts- 
In order to correctly incorporate all effects at next-to-leading order in the chiral expansion, it is necessary, for those 
terms linear in ir 3 , ir 8 and not proportional to the fourth order LEC's £6,7,8 m the representations of the currents, to 
use a transformation between the tt 3 , tts and tt°, r\ bases which is valid to 1-loop order. The relevant expressions are 
[16]: 



_ 7 -l/2 



7T = Z. 



n 



7T3 + #l7T 8 



-1/2 ^ 



(15) 



where Z 3 and Z% are the 7r 3 and tt 8 wavefunction renormalization constants to 1-loop order in the isospin symmetry 
limit and 



1 7 4 

1 - T^Mtt + + Mr; 



3to 2 + m 2 
64tt 2 / 2 



mt 



■ \og(ml/m 2 K ) 



2 2 



#2 — &o 



64tt 2 / 2 
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- ) (l + log(m% /m 2 )) - 32( ^ 2 ml) {ZU 7 + LI) 



f 2 



1 - y Mtt + g m + Vr, + 



3m 2 



64tt 2 / 2 



1 + 



m 2 K — m 2 . 



\og{m 2 Jm 2 K ) 



2 2 

— m 



64tt 2 / : 



i ] (1 + log^/M 2 )) - 32{7h \ 2 ml \ 3LL , LI . 



(16) 



In Eqs. (16), fip = [m 2 p log(m 2 P / /j, 2 )]/32it 2 f 2 , with [i the (dimensional regularization) renormalization scale, L r k {n) are 
the renormalized version of the fourth order LEC's, and m 2 K is the average of the neutral and charged kaon squared 
masses. Standard manipulations then lead to the following expressions for the vacuum-to- ir°, rj matrix elements, valid 
to 0(md — m u ), and to fourth order in the chiral expansion: 



< 0\P U - P d \ir° >= 



rh 

,2 



<0\P u + P d \n° >=^^ 
V3m 



1 + 



mZ 



m 2 K -ml 



log 



in 



K 



f 

< 0\P,\n° >= - 



288i?o / ~ \ T r 
——(m s - m)L 7 



+ 



\/3m 

B a (m s + m) 
16tt 2 / 2 



-7^ + Ari K + 'ijJ-r, 
P 

16tt 2 / 2 { + 



—^-{m s - m)L s 



+ 2\i K \ B (m s - m) 



i 3m 2 + ml 
y 64tt 2 / 2 

log^/M 2 )) 



f 2 



16tt 2 / 2 



(1 + \og(m 2 K /^ 2 )) 



{m\ - ml) 



log -f 



in 



K 



(17) 



(18) 



(19) 
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<0\P u -Pd\v> 




-7fJ, n + 6fJ, K + Hr) 



3/ 2 

,2 „_2 \ 



3m 2 v + ml\ f 



( m n ~ ml) , 
64tt 2 / 2 ' 



I 64tt 2 / 2 
(l + log^/M 2 )) 



< 0|P„ + P d |r, > = 



\/3m 



96P , Mrr 32S 
— (m s - m)L 7 ■ 

2/^ + 2/^A 16B 



3 p-( m s ~rh)L r 5 —(m s - m)L r 7 ^-(m s - m)Lg 



/ 2 



(20) 

(21) 



< 0|P s |ry >= 



V3m 



i - ( — ^ — 1 - -jjr( ms - m > L 5 + —p-y ms - m > L 7 



H — j2-(m s - m)L 8 



(22) 



The cancellation of divergences, and the overall scale-independence of the expressions for the matrix elements above, 
provide non-trivial checks on the calculations, as do the relations 



< Q^A^P >=< 0\2m,fPf\P > + anomalous terms 



(23) 



which follow from the Ward identities between the divergence of the axial currents and the pseudoscalar currents. 
(Note that the isospin-breaking vacuum-to- ir° , r/ matrix elements of the axial currents Ag and A% have been evaluated 
previously by Gasser and Leutwyler [26], and that the anomalous terms in Eq. (23) cancel when one considers the 
combinations which enter A^ and Ag. The vacuum-to-pseudoscalar matrix elements of the anomalous terms are given 
explicitly in Ref. [29].) In bringing the expressions for the matrix elements into the forms displayed above, we have 
made use of the 1-loop expressions for m 2 and [26]. The appearance of the common multiplicative factor m 2 / w /m 
in all the matrix elements has no physical significance; the equations are simply written in this form in order to 
simplify the later calculation of the mixing parameters relevant to the "pure" n° interpolating field. 



III. INTERPOLATING FIELDS FOR THE NEUTRAL PION 

As explained above, the MH tt° interpolating field choice Eq. (3) is not unique. In fact, if we consider 

P{6) = (P u - P d ) + e{6) (cos9 P 8 + sine P ) (24) 

where P = (P u + P<i + P s )/V3, the results of Eqs. (18), (19) show that, in general, 

< 0|cos(9P 8 +sin6»P |77 >^ . (25) 

As a result, one can, using the construction above, chose e(9) in Eq. (24) in such a way that < 0|-P(#)|?? >= for any 
value of e. Any such P(0) can then be employed in an analogue of the MH analysis of the three point correlator, in 
the sense that, in all such analyses, the unwanted r\ contributions will have been removed by the corresponding choice 
of e(0). 

The existence of the family of 7r° interpolating fields above is potentially useful as a means of investigating features 
of the three-point-function method that have been the subject of criticism in the literature. The most serious such 
criticism of the method as a way of estimating the ttNN couplings is the objection that it requires making the 
assumption of pion dominance in the pseudoscalar channel all the way out to values of q 2 of order — 1 GeV 2 for which 
one can begin to hope that the operator product expansion (OPE) analysis of the correlator might be valid [8]. The 
danger, of course, is that, so far from the pion pole, the contributions of other resonances (77', 77(1295), 7r(1300), etc.) 
might no longer be negligible. Although MH have not attempted to investigate this aspect of the method in detail, 
they argue that their choice of ir° interpolating field (corresponding to e — 0, and the leading chiral order estimate 
for e(6) above) represents the unique one that, at least to leading order, also has no overlap with the 77'. At first sight 
this seems a plausible argument since, in the SU(3)f and isospin symmetry limits, P u — Pd couples only to I = 1 
states and P 8 only to flavor octet states. On closer inspection, however, one can see that < 0|P(0)|r/ >^ [30], 
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so that there remain potentially dangerous higher resonance contaminations. Although this point is part of a more 
detailed analysis of higher resonance contaminations in the three-point-function approach, the results of which will 
be presented elsewhere, we will briefly outline the key aspects of the argument below. The main point of relevance to 
the present discussion is that, although one can, in principle, find a value for 8 for which at least the vacuum-to-// 
matrix element vanishes, the information required to actually determine this value is not available. Thus, for any 
particular choice of 8, one has a coupling of P{8) to the rf of unknown strength (in addition to potential couplings to 
other higher pseudoscalar resonances), which presents a serious problem for the three-point-function method. 

One way of dealing with this problem is to investigate the extracted values of the isospin-breaking nNN couplings 
as a function of 8. If one finds a region of 8 in which the results are not sensitive to variations of 8, then one will have 
some a posteriori justification for assuming that the channel of interest is, indeed, dominated by the lowest-lying (ir°) 
state. Note that the coupling of P(8) to the physical 7r°, from Eqs. (17)— (19) above, is, to C(md — mj, m 2 / 2 /m, 
independent of 8. The vacuum-to- 7r° , rj matrix elements above provide the essential input to such an analysis, and we 
will recast them in a form suitable for use in such an analysis below. 

Before doing so, however, let us briefly outline why the MH field choice, 8 = 0, does not actually remove the 
coupling to the r/. Recall that the quark mass matrix responsible for flavor and isospin breaking can be decomposed 
as 



M 



1, 
- 2 {m d - 



m u )X 3 - —^( m s 



(26) 



To leading order in the isospin-breaking and flavor-breaking mass differences, therefore, the vacuum-to-?/ matrix 
elements of P u — Pd and P§ are related simply by the ratio of mass-dependent coefficients of A3 and Ag in Eq. (26). 
It then follows straightforwardly [30] that, to leading order in the masses, 



< 0|P(0)|r/ >= 30 o < 0|P 8 |f/ > 



(27) 



where < Ol-Psl 7 / >= 0(m s — rh). The RHS is thus non-zero, and 0(rrid — m u ). If one considers 8 7^ 0, of course, one 
has, in addition to the 8f x 8f — > If reduced matrix element which governs < 0|P U — Pd\v' > an d < 0|Ps|?/ >> the 
If x If — > If reduced matrix element relevant to < 0|Po|r/ >• One can then certainly, in principle, find a 9 such 
that < 0\P(8)\rj' >= 0. The problem is that, even to do so at leading order in the quark masses, one would need to 
know the ratio of the two reduced matrix elements above, and this information is not available. Moreover, even if it 
were, this would not necessarily ensure that, for such a value of 9, the couplings of the higher resonances were small 
for the same value of 8. We will return to these issues, and to a discussion of the reliability of the three-point-function 
method, in more detail in a future publication [30]. 

Let us return, then, to our 1-loop evaluation of e(8). From Eqs. (7), (20), (21) and (22) it follows that 



e(8) 



cos 8 



N 



(28) 



where 



N = 



/ -7^ + 6/ig+^A ( 3w 2 + ml \ / 
I 3/ 2 J + [ 64tt 2 / 2 ) \ 



"K 7T 



log 



D 



(to 2 - ml) IQB 0( 
+ , ; o ,0 (1 + logCmJf/V)) - - 177T (m s - m)Ll 



96B 



3.f 2 



- r 32Pn 



K 



(m s -m)L r 7 —{m s -m)Ll 

'it&nd ^ + 2(1 + tan6>)^F; + (tan# — 2)fi n \ l6B (m s — rh) 

1 + I TTl^ I Lie. 



3/ 2 



32B Q (m s -m) 32B (m s - rh) 

■tan#L 7 —2 (tant* - 2)L 8 



3./ 2 



P 



3/ 2 



(29) 



(30) 



The factor 8 / cos8 represents the leading order contribution to e(0) in the chiral expansion and, of course, reproduces 
e MH (corresponding to 8 = 0) at this order. The remaining factors represent the next-to-leading order corrections. We 
will discuss their magnitude below. One should note that, as 8 — > ir/2, 3tt/2, the denominator of Eq. (28) tends to a 
finite value involving the chiral log and fourth order LEC terms. This means that, near such values of 8, e(8) becomes 
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large, and with significant uncertainties (associated with uncalculated next-to-next-to-leading order contributions to 
the denominator). As such, the corresponding interpolating field choices, though in principle reasonable, are not 
useful practically. In order to avoid these regions of 9 it is convenient to consider a re-parametrization of the family 
of interpolating fields in which one writes 



P(a) = (P u - P d ) + e(a) a(P u + P d ) + (1 + a)P t 

1 



{Pu - Pd) + e(a) 



V3 



-P 8 + (3a + 1)P 



where, as usual, e(a) is to be chosen for each a in such a way that 

< 0\P(a)\ri >= . 



(31) 



(32) 



The second line of Eq. (31) shows that such a form for P(a) represents a simple reparametrization of P(0), in which 
a factor of — \/3cos# has been absorbed into e and tan 9 has been rewritten as 3a + 1. The regions of 9 we wish to 
avoid then correspond to a — > oo. One then easily shows that Eq. (32) implies 



e(a) = -V3<9 



1 

32 



-(10 + 9a) + 6(1 + a)fi K + (4 + 3a)^ 
3/ 2 



3/2 -(4 + 3a) (ml - ml) (3L? + LJ) + (^^J (l 
-<S)) + ^( 1 + 1 o g (^)) 



= - V36» [1.299 + 0.0894 a] , 



(33) 



where, in arriving at Eq. (33), we have made an expansion of the denominator in Eq. (7) so that the whole expression 
in now explicitly consistent to 1-loop order. The reparametrization has been chosen such that, at leading order in the 
chiral expansion, e(a) reduces to — V3 9 0} so the leading 1 in the square brackets in Eq. (33) represents the tree-level 
contribution. The remaining terms are the 1-loop corrections and, as expected, are not small. This can be seen from 
the last line of Eq. (33), where the deviation of the term 1.299 + 0.894 a from 1 given the size of the 1-loop corrections. 
For reference sake we note that the MH choice 9 — corresponds to a = —1/3, and the choice in which there is no P s 
content to the interpolating field corresponds to a = — 1. (To completely convert back to 9 notation for the MH case, 
one must also absorb the constant — 1/\/3 back into the definition of e.) The enhancements due to next-to-leading 
order corrections are, in these cases, 1.27 and 1.21, respectively. Amusingly, the choice a = —4/3 removes the LEC 
contributions completely from the expression for e. 



IV. SUMMARY 



We have evaluated the leading 0(m,d — m u ) strong isospin-breaking contributions to the matrix elements < 
0\Pf\ir°,r) > at next-to-leading order in the chiral expansion. The results are given in Eqs. (17)-(22). We have 
also showed how to use this information to construct a family of ir° interpolating fields which could, in future, be 
employed to investigate the reliability of analyses which attempt to extract the isospin-breaking ttNN couplings via 
the three point function method in QCD sum rules. The isospin-breaking parameters, e(a), required in this construc- 
tion in order to remove unwanted r\ contaminations from the extracted tt° couplings, were shown to receive significant 
contributions at next-to-leading order. As a result, if one wishes to remove such contamination, it is essential that 
the full next-to- leading order expressions, given in Eq. (33), be used. 
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FIG. 1. Graphs contributing to the vacuum-to-7r , 77 matrix elements. The wavy lines represent the pseudoscalar currents, 
the solid dots the vertices generated by the 0(q 2 ) part £ e //, the 'X' and the circle enclosing an 'x' vertices generated by the 
0(q 4 ) part of £ e ff- 
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